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We consider the effect of modified gravity on the peak structure of the Cosmic Microwave Back- 
ground (CMB) spectrum. We focus on simple models of modified gravity mediated by a massive 
scalar field coupled to both baryons and cold dark matter. This captures the features of chameleon, 
symmetron, dilaton and f(R) models. We find that the CMB peaks can be affected in three inde- 
pendent ways provided the Compton radius of the massive scalar is not far-off the sound horizon 
at last scattering. When the coupling of the massive scalar to Cold Dark Matter (CDM) is large, 
the anomalous growth of the CDM perturbation inside the Compton radius induces a change in the 
peak amplitudes. When the coupling to baryons is moderately large, the speed of sound is modified 
and the peaks shifted to higher momenta. Finally when both couplings are non-vanishing, a new 
O , contribution proportional to the Newton potential appears in the Sachs- Wolfe temperature and in- 

• creases the peak amplitudes. We also show how, given any temporal evolution of the scalar field 

mass, one can engineer a corresponding modified gravity model of the chameleon type. This opens 
up the possibility of having independent constraints on modified gravity from the CMB peaks and 
large scale structures at low redshifts. 
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■ Modified gravity is a tantalising idea which has received a lot of attention recently, sec [1, 2] and references therein. 
| This is partially due to the quest for an explanation to the acceleration of the Universe. Another reason is the 

l— ~~ 1 ■ possibility of constructing meaningful models of massive gravity. In most of these attempts to modify gravity, scalar 
| fields play a prominent role. From the dark energy point of view, scalar fields are plagued with the fifth force problem 
^ . whereby the scalars generating the acceleration of the Universe would lead to an extraneous contribution to Newton's 
' law. This can be cured by screening mechanisms locally. Indeed, dense regions such as the sun or the earth can trap 
\Q , the scalar field and a minute deviation to Newton's law ensues. This the case of the chameleon models [3-7]. In other 
OO ■ cases such as the galileon [10], symmetron [8, 11] and dilaton [9], the coupling to matter almost vanishes in dense 
environments. 

Q\ ■ Unfortunately, these models behave like A-CDM at the background level at least since Big Bang Nucleosynthesis 
jeopardising any attempt to distinguish them from the cosmological standard model. Fortunately, they modify gravity 
,— 1 . at the perturbation level and lead to observable effects on the growth of structures in the recent past of the Universe. 

Typically, gravity is modified inside the Compton radius of the scalar field which must be of order of a few Mpc now. 
J> , As the growth of structures is observed using galaxies as tracers of Newton's potential, the analysis is hampered by 
K> ' ^he galaxy bias and therefore subject to some uncertainty. On the other hand, the physics of the CMB is linear and 
m | evades all non-linear problems. The CMB could serve as a clean template for the seach of modified gravity effects. 

■ This is what we will pursue in this paper. 

Modified gravity will be most effective at last scattering when the Compton radius is close to the sound horizon. 
In this case, scales inside the sound horizon would be affected by a modification of gravity. An extreme case, where 
the mass vanishes and all scales are affected by modified gravity, has been considered in the Brans-Dicke framework 
in [12, 13] and recently revisited by [14] where it has been found that the constraints are less severe than the Cassini 
bound on the Eddington parameter [15]. Here we will always consider the massive case with a Compton radius in the 
range of the CMB peak scales. In this case, three effects can occur. The first one is driven by the coupling to CDM 
and consists in an anomalous growth of Newton's potential. This increases the peak amplitude. Another effect results 
from a decrease of the speed of sound coming from the coupling to baryons. If this coupling is too large, an instability 
with an imaginary speed of sound sets in as noticed in [16]. Finally a more subtle effect comes from both couplings 
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and induces a modification of the contribution of Newton's potential to the Sachs- Wolfe effect [18, 19]. 

We also show that the knowledge of the scalar field mass around the last scattering surface is enough to reconstruct 
some part of the interaction potential of a chameleon model with such a mass and couplings. If some knowledge of 
the mass function at later times could also be inferred from large scale structures, a larger range of the interaction 
potential could be rebuilt. Eventually, if the mass function could also be extracted from data in the intermediate 
regime corresponding to the cosmic dark ages, an almost complete picture of the interaction potential of the scalar 
field could be drawn. Together with the values of the couplings to baryons and CDM which could be obtained from 
the deviations of the CMB peaks from their A-CDM values, an almost complete description of the source of modified 
gravity since the last scattering epoch could be obtained. 

In the first part, we analyse linear perturbations and the growth of structures in a modified gravity context. Then 
we use the tight binding approximation [21] to study the effects of modified gravity on the CMB peaks. Then, we 
show that one can always engineer any modification of gravity on large scales by choosing the time evolution of the 
scalar field mass. Finally we discuss our conclusions. 

II. MODIFIED GRAVITY AND LINEAR PERTURBATIONS 
A. Perturbations 

Modified gravity will be modelled using a massive scalar field <p with a time varying mass m(a) and coupling 
constants /3& and j3 c to both baryons and Cold Dark Matter (CDM). We will not assume that the couplings are 
constant, as some interesting models such as the environmentally dependent dilaton have field dependent couplings 
[9]. On large cosmological scales and at the perturbative level, models of modified gravity such as f(R) theories, 
chameleons, dilatons, symmetrons and galileons all involve a scalar field with couplings to matter [4, 7-11]. On 
shorter distances, one must invoke the existence of a screening mechanism in order to evade tight gravitational tests 
in the solar system, such as either the chameleon or the Vainshtcin screening properties. At the level of the CMB, 
the screening of gravity on small scales is not present and large effects may ensue, depending on the range of the fifth 
force mediated by the scalar field and the intensity of the coupling to matter. 

We focus on scalar fields in the Einstein frame where the Einstein equations are preserved and depend on the energy 
momentum tensor of the scalar field: 

T* u = d^d u <t>-g^(dct>) 2 +V). (1) 
Matter couples to both gravity and the scalar field via the metrics 

ajiu = A l{4>)g»u, i = b,c (2) 

for each matter species. Their couplings to matter are defined as 

= TO P1 — (3) 

which may be field dependent in dilatonic models for instance. It is a universal constant j3 — ^= in f(R) models. 

In all the models that we will consider, the background cosmology is tantamount to a A-CDM model since Big 
Bang Nucleosynthesis (BBN). Deviations from General Relativity only appear at the perturbative level. Moreover, 
the effective potential for the scalar field in the matter era is modified by the presence of matter as a consequence of 
the non-trivial matter couplings 

v eS {4>) = v{4>) + Y, A i{4>)Pi (4) 

i 

where the sum is taken over the non-relativistic species and pi is the conserved energy density of the ith fluid. The 
potential acquires a slowly varying minimum <fi(pi) which is an attractor as long as the mass m 2 p = U(p<) ^ s 

larger than the Hubble rate [5]. We will always assume that this is the case in the following. Due to the interaction 
with the scalar field, matter is not conserved and satisfies 



(5) 
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where T? v is the energy momentum tensor of the ith species, T its trace and k 2 = 8itGn- Notice that radiation is not 
affected by the presence of the coupled scalar field. The Einstein equation is not modified and reads 

- \g»uR = 8ttG w (^ T% + T* v ) (6) 

i 

where rrip 2 = 8ttGn- The conservation equations and the Einstein equations arc enough to characterise the time 
evolution of the fields. 

We are interested in the first order perturbations of the Einstein equation and the conservation equations. In the 
absence of anisotropic stress, the metric can be described in the conformal Newton gauge 

ds 2 = a(r/)(-(l + 2<S>)drf + (1 - 2<P)dx 2 ) (7) 

where $ is Newton's potential. In this gauge the relevant Einstein equations involve 5T® and 5Tq . For the scalar field 
we find that 

<5T° = (8) 

where ' = d/drj and 

ST g = ^$ - ^ - V+St (9) 
which implies that the Poisson equation becomes (in Fourier modes) 

fc 2 $ + 3H(<P' + H&) = -AnG N a 2 V So, - 47rG A ra 2 (^ T $ - ^5d>' - VaS6) (10) 

and 

fc 2 ($' + U<&) = AnG N a 2 ^(1 + Wi)0i + AirG N k 2 4>'8(t> (11) 

i 

where the Hubble rate is H = a' /a and the equation of state of each species is w;. We have defined the divergence of 
the velocity field of each species 9 = diV % . The perturbed Klein-Gordon equation can be expressed as 

5<j)" + mH' + {k 2 + a 2 m 2 )5<t) + 2$a 2 V off0 - 4$>' + K A a 2 {p c A c 8p c + f3 b A b 8p b ) = (12) 

where Spi = piSi. We have denoted by p b the baryon energy density and p c the CDM energy density. 
The conservation equation for the CDM particles lead to the coupled equations 

5' c = -(9 C - W) + K^<t>'5<t> + n 4 p c A c S<j)' (13) 

and 

9' c = -Wo + fc 2 $ + mk 2 f3 c A c 6<t> - p c A c 4>'e c . (14) 
The baryons are coupled to the photons by Thompson scattering implying that the conservation equations are modified 

S' b = -9b + 3$' + K 4 fy<f>'6(j) + f3 b K 4 A b 5(j>' (15) 

and 

9' b = -H9 b + ^^(0 7 - b ) + fc 2 $ + p b A b K A k 2 8dp - pM9 b . (16) 

The Thompson scattering cross section depends on the fine structure constant a and the electron mass m e . The 
electron mass is scalar field dependent due to the conformal rescaling of the metric, i.e. m e is proportional to Ab(<f>). 
Perturbations of the scalar field would have an effect at the level of second order perturbation theory. At the linear 
level, the dependence of the electron mass on must only be taken into account at the background level. As the scalar 
field tracks the minimum of the effective potential, the time variation of the mass is electron mass in one Hubble time 
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is suppressed by O(^) <C 1, see (24), and can therefore be neglected. If the fine structure constant depends on the 
scalar field, its time dependence can also be neglected for the same reason. 
We will identify the speed of sound in the absence of modified gravity as 



9 1 1 

c ° = 3 IT* (17) 



where the baryon to photon ratio is 



< 18 > 

which is about 0.6 at last scattering. Finally we shall work in the fluid approximation for the photons. The Boltzmann 
hierarchy is not altered by the presence of a scalar field. We find 

5', f = -1ft, (19) 

and 



^ _ X ,1.2 



— <5 7 + r$ + an e o- T (6'&-6» 7 ). (20) 



This completes our description of the perturbation equations. 

We also need to specify the initial conditions for all the perturbations. We will be interested in modes which will 
enter the horizon before radiation-matter equality. Adiabatic initial conditions are determined by 



This is related to the initial Newton potential as 



% = S° b = (21) 



S° b = -|*o- (22) 



We will always express the perturbed quantities in terms of the initial Newton potential. 



B. Growth of structures 



The growth of CDM perturbations is affected by the modification of gravity due to the scalar field. Let us consider 
sub- horizon perturbations where k 3> H- We assume that the background configuration of the scalar field tracks the 
minimum of the effective potential: 

= -n A l3 c A c p c - K 4 /3 b A b p b . (23) 
This implies that we can also neglect terms in <fr' as 

^f= 9 J2A^~ (24) 

i 

where m 3> H. In the sub-horizon limit, we also neglect the time variation of Newton's potential as terms in /c 2 <3? 
dominate over $. In this case the Klein-Gordon equation becomes an algebraic equation 

(i c A c 8p c + /3 b A b Sp b 

0(p = -K4 p • (25) 

K +m 2 

As /3 b can differ from /3 C , the baryon contribution to the scalar field perturbation cannot be neglected. Similarly the 
Poisson equation becomes 

fc 2 $ = -47rG A ra 2 ((5p c + ( 5 /9f) + Ap 7 ). (26) 
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We can neglect the contribution of the baryon and photon perturbations 8b 1 = 6f>b '' 1 to Newton's potential in the 

matter era as the CDM density contrast S c = grows rapidly. Similarly in the radiation era, the baryons and 
photons are tightly coupled and their oscillatory perturbations average to zero. As a consequence we can neglect the 
baryons and photons in the Poisson equation up to the last scattering time. The CDM equations become 



and 



leading to the growth equation 



S'c = -0c (27) 
-n9 c + k 2 ^ + f3 c A c k 2 K4,S(j) (28) 



5'; + H5> c -Ih 2 ^ (1+ 2p2cA } 2 ) = i) (29) 



When A c ~ 1 and p c 3> p&,/9 7 , this is an equation which can be easily integrated deep inside the Compton radius 
k/a 3> m and outside the Compton radius fc/a < ra in the matter dominated era. Outside the Compton radius, the 
CDM density contrast grows like: 

5 C ~ a (30) 

as in General Relativity while inside the Compton radius the growth of structures is enhanced due to the factor 
(1 + 2/3;?) corresponding to an increase of Newton's constant 



-l + -i/l + 24(l + 2/3 2 ) 

6 C ~ a 1 (31) 

In particular this implies that Newton's potential is not constant in time on small enough scales. A well known 
consequence of this fact is a new contribution to the Integrated Sachs Wolf (ISW) effect for the CMB spectrum, 
see [18] and references therein. Finally, during the radiation era the baryonic and photonic perturbations are tightly 
coupled and oscillates while the CDM perturbations grow logarithmically as 

5 c = 5 (l + 2^^1n-) (32) 

where and Hk are the time when the momentum k enters the horizon, when Hk = k, and fc cq is such that scales 
with k > k eq enter the horizon before matter-radiation equality. It turns out that jr'-ak = o, eq implying that 

<S c = <5 (l + 21n— ) (33) 

where <5o is the density contrast at horizon entry. Initially we have that So = — 1$(0) where $(0) is the constant 
Newton potential outside the horizon and we have taken <E>| afc =0 which implies that \ a k = 2. 

The linear approximation can be used provided the screening of the scalar field effects is not present. The thin-shell 
mechanism occurs when the variation of the scalar field profile due to an overdensity satisfies 

Ki5(j) < 6/3 c $ (34) 

where we simplify the problem and assume that only CDM is responsible for the overdensity. Using (25) and (26), 
this is equivalent to 

k <C am{a). (35) 

In the following we will be interested in scales which are close enough to the Compton wavelength, implying that the 
linear analysis is valid. 

At the beginning of the matter era, the growth of perturbations depends on the behaviour of am(a). When am(a) 
increases, the scale k can already be inside the Compton radius at matter equality and then grow anomalously before 
leaving the Compton radius. Schematically, this leads to a growth in In a, followed by a"/ 2 , followed by a. If the scale 
k is not in the Compton radius at equality and am(a) increases, it will never be. On the other hand when am(a) 
decreases, if the scale k is not in the Compton radius at equality, it will enter the Compton radius at a later time 
implying a growth in In a, followed by a, and then a u / 2 . These patterns will be explicit in the examples we will present 
in the next section. 
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C. From horizon entry to last scattering 

We will give an explicit example of anomalous growth of structures from horizon entry to the last scattering surface. 
This will make explicit the general remarks of the previous section and will be useful when we discuss the CMB in 
the tight coupling approximation [21]. We focus on models where the time dependence of the mass is given by 

m(o)=mi(— ) r (36) 
aLS 

where aLS is the scale factor at last scattering and r is an index. The scale mi will be chosen to be of the order of the 
(inverse) CMB length scales, i.e. close to the sound horizon (see sections III-C and III-D for a description of models 
leading to this power law behaviour ). The comoving Compton length increases when r > — 1 and decreases when 
r < — 1. This implies that scales enter the Compton radius when r < — 1, and they eventually leave the Compton 
radius when r > — 1. More precisely, scales coincide with the Compton radius when 

acompton = a L s(^) 1/{r+1) (37) 

where fci = aigmi. When r < — 1, scales k < k\ enter the Compton radius after last scattering while scales k > ki 
enter the Compton length before last scattering. When r > — 1, scales with k > k\ leave the Compton radius after 
last scattering while scales k < k\ leave the Compton radius before last scattering. 

Let us consider a given scale k which enters the horizon before matter-radiation equality. Up to matter-radiation 
equality, the density contrast 5 C grows logarithmically. After the equivalence, one must distinguish whether am (a) 
increases or decreases. 

1. am(a) decreases 

This corresponds to r < — 1. It is very useful to determine when scales enter the Compton horizon after matter- 
radiation equality This happens when 

k < k c = fcl (i±^2-)-(r+D (38) 

1 + Z L S 

Hence when k eq < k < ki, scales enter the Compton length after last scattering. Here k eq is such that modes k > k eq 
enter the horizon before matter-radiation equality. Such a mode grows logarithmically till matter-radiation equality 
and then undergoes a normal growth implying that 

%£ = ^£ ( l + 2 l n ^£) (39 ) 
d a eq a k 

When ki < k < k c , a mode first grows logarithmically till matter-radiation equality, then grows at a normal rate 
when outside the Compton radius and then grows anomalously. The density contrast at Compton radius entry is 

|c = acompton (1 + 21n fle i) _ (4Q) 

oo a eq a k 

The growth is anomalous between acompton and clls implying that 

6ls_ = £W(A)(i-"/2)/<r+i)( 1 + 2 in£2i). (41) 
Sq a eq fci ak 

Notice the anomalous power law dependence. When k > k c , the density contrast is simply 



5ls 
So 

which corresponds to an anomalously enhanced density contrast compared to General Relativity 



^S =( ^ r /2 (1 + 21n ^ i) (42) 
t>0 a eq &k 



7 



2. am(a) increases 

This corresponds to r > — 1. Scales leave the Compton horizon after matter-radiation equality when 

k > k c = ki( \ + ZLS ) r+1 (43) 

Hence when k eq < k < k c , a mode first grows logarithmically till matter-radiation equality and then grows at a normal 
rate until last scattering as it is outside the Compton radius. On the other hand, when k\ > k > k c modes grow 
logarithmically, then anomalously and eventually at a normal rate when they have left the Compton radius before 
last scattering. The density contrast when leaving the Compton length is 

|c = ( aCompton ^ /2(1 + 21n ae i) _ (44) 

do a eq at 

The growth is normal between acompton and a^s implying that 

S_LS = OLS ( _fc )W2 -l)/(r+l) (1 +21n ^i). (45) 

do a eq ki at 
Notice the anomalous power law dependence. Finally when k > k±, scales grow logarithmically and then anomalously 

%£ = ( ^£ r / 2( l + 21n^) (46) 
do a eq afc 

which corresponds to an anomalously enhanced density contrast compared to General Relativity. 

3. Discussion 

We have deduced that the growth of the CDM density contrast is affected in an interval due to the fact that scales 
enter/leave the Compton radius. When such an anomalous power law is present, we find therefore that Newton's 
potential is not growing logarithmically only but 

fc 2 $oc( fc ) W2-l)/(r+l) (1+21n ^e i) (4?) 
fcl <2fc 

which is due to the modification of gravity in the [fci,fc c ] (respectively [fc c ,fci]) interval. In the numerical examples, 
we will take r = 1 and r = —2. 

III. THE COSMIC MICROWAVE BACKGROUND 

A. The tight binding approximation 

We are interested in the temperature fluctuations of the CMB as given by the Sachs- Wolfe formula 

ST 1 f° 

— = -5 7 + $ + e iV \ + 2 J e <t>'d V (48) 

where is a vector along the line of sight. Effects of modified gravity are multiple. As we have already mentioned, 
the fact that Newton's potential varies in time leads to a contribution to the ISW effect, see for instance [18] and 
references therein. Here we arc interested in intermediate scales where the peak structure of the CMB is relevant. We 
will study the effect of modified gravity on the temperature anisotropy 



6 = -<5 7 + $ (49) 



using the tight binding approximation [21], allowing us to capture the essence of the consequences of modified gravity 
on the CMB peaks. This approach is not precise enough to impose tight bounds on the couplings. We will only use it 
to analyse the possible consequences of modified gravity on the CMB qualitatively. 
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On sub-horizon scales and neglecting the time variation of <f> we have 

s' b = -e b (50) 

and 

S!y = (51) 

In the tight binding approximation, the photon and baryon density contrasts are linked by 

an e a T {e b -e i ) = 9' 1 - —8 1 ~k 2 <S> (52) 

which does not involve the scalar field directly. When an e (jT is larger than the Hubble rate this implies that 

6 b ps 7 (53) 

and therefore 

S b « \5, (54) 

5 b + TT-5 S 'b + c ^ = - kH ~ ^PbA b k 2 -^—5(l>. (55) 
1 + ti H + 1 

It is particularly useful to define 

S b = (l + R)- 1 / 2 S (56) 



leading to 



from which we deduce that 



5" + c 2 k 2 S = -k 2 (l + R) 1 ' 2 ® - K4 f3 b A b k 2 —JL— 5cl> (57) 



in the sub-horizon limit. Neglecting the baryons and photons in the Poisson equation we find that 

U& R 
i£ R+l 



S" + c 2 k 2 (l - ^f Al ™ 2 )5 = -k 2 (l + Rf' 2 {l + 2 ? bA ±j b ~ - )*. (58) 



Finally we can introduce 



S^^ + ny^SJ^^* (59, 



where the effective speed of sound is 

r 2 =^ 2 n - - 

k 2 + m 2 a 



e.^-m***.) ( 60) 



which leads to 

5" + c 2 J = 0. (61) 

The effective speed of sound can become negative when (3 b is too large leading to an instability already noted in [16]. 
The WKB solutions of this equations are 

5 = c~ 1/2 (Asmkr s + Bcoskf s ) (62) 

where r s (rf) = c s drj is the modified sound horizon. We have taken the initial condition to be at r]k when scales 

becomes sub-horizon. % = k. The WKB approximation is valid as long as w' <C w 2 where w = kc s . This is valid 
inside the horizon. Initially both 8 1 and <E> are constant implying that 

(5 = 3 1/4 (l + i?) 1/4 Bcosf 5 Jc (63) 
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and we find that 

B= (1+ 3 f / f / \ l + 2R k ) 3 -m (64) 

where we have set the initial conditions as adiabatic at horizon entry in the radiation era. This leads to 

5= ^(l + i? fc ) 1/4 (l + 2i? fc )*(0)(l + i?) 1/4 cosr s fc (65) 

where we have used Sb{0) = — 1$(0) outside the horizon. We have also used the fact that initially all scales are outside 
the horizon and thus outside the Compton radius. It is important to relate S to 

9 = 3(1 + R)^ 3c ' + 1 + R + l ] W (66) 



This can also be written 



where we have defined the effective couplings 



ho = (68) 



fc 2 



There are three sources of alteration to the peak structure. First of all the speed of sound is modified due to /?{,. Then 
the growth of Newton's potential is also different from the one in General Relativity due to (3 C . Finally, the amplitude 
of the term in $ appearing in the expression for is also modified with a contribution depending on /3b/3 c . 

B. Phenomenology 

Assuming that the equation of state varies abruptly at equality we can write 

o = a eq ( o ) + 3v \/ ] -^ + l) 2 /^ 3 ™) (69) 
eqy 3(l + w) {l + z eq )t eq ' v 1 



where r\ is the conformal time. Defining 



2 I 

k eq (w) = -— r-. r — (70) 

eq v ; 3(1 + w) (l+z eq )t eq v ' 



in both the matter and radiation eras, we find that scales enter the horizon when Hk = k at 

°fc _ f k eq (w) ^ 2 /3(l+w) 

a eq h 



(71) 



corresponding to 



2 1 1 

Vk = veg + -t i o ( y _ i — 7-5) ■ ( 72 ) 

1 + 3w k K eq (w) 



Vk = Veg + (-- T-^) (73) 



Scales with k > k eq (w) enter the horizon before equality at 

,1 1 
'A: A; e9 (w) 

We consider scales entering the horizon in the radiation era and put k eq = fc eg (l/3), then we have 

2(1+ Zeq)t eq 
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where 



" (1 + X)3/ 2 



Scales corresponding to the CMB peaks arc given by 



where a good approximation is 



or equivalently 



and 



r s {VLS) 



r s (r)Ls) ~ ^(vls - Vk) (77) 



r S (VLs) ~ 7^{VLS - Veq + ^ - -£-)) 



<■-</ 



ms-iM^j-ii^-^) 1 ' 2 -!). (79) 
k eq 1 + z LS 

The corresponding Compton wavelength is 

m p = (1 + z LS )k p . (80) 
It is convenient to compare this to the Hubble rate (neglecting the effect of the late acceleration of the universe) 

H LS = (1 + z LS ) 3/2 H (81) 



^-(I + zls)-^. (82) 



tt2 ^ 1 rr2 

W L5 ""0 



Numerically we find with these approximations that ki « 0.035 (h x .Mpc) 1 and 



,2 



H LS 



12.8p z . (83) 



This is not satisfied for dilatons as m 2 / H 2 = 0(A 2 ) is very large [9]. As f3 is small for dilatons, there should be no 
CMB effect in this class of models. 

In the following, we set the initial conditions for a mode k at rjk- In the figures, we have selected wavenumbers 
corresponding to the first two peaks and chosen a power law dependence for the mass of the scalar field 

m{a) = m c ( — Y (84) 

where m c is a characteristic scale depending on the model. With this parametrisation, the mass is always greater than 
the Hubble rate since equality provided 

r + g < " L z 3 . (85) 
2 ln^i 

l+ZLS 

Scales outside the Compton radius enter the Compton radius when 

k = a k m(a k ) (86) 

after horizon entry. This implies that for a given k, gravity is modified as soon as m(a)a < k. We find that the 
Newtonian potential grows logarithmically in the radiation era leading to a small growth of fc 2 $ when no modification 
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FIG. 1: Newton's potential fc 2 $ at last scattering for /3b — 0, /3 C = 100, r = —2 (left) and r = 1 (right) as a function of fc in 
(h _1 .Mpc) _1 compared to the case with no modification of gravity. The logarithmic increase is largely modified due to the large 
value of /3 C . Here we have fci ~ 0.035(h _1 .Mpc) _1 , and for r- 2, k c ~ 0.1(h _1 .Mpc) _1 , while for r = 1, k c ~ 0.006(h _1 .Mpc) -1 . 
In both cases, the anomalous growth of perturbations due to modified gravity can be seen for scales between ki and k c , 
respectively k c and ki. For smaller scales, the growth is only logarithmic. 
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FIG. 2: Newton's potential fc 2 <3? at last scattering for f3 — 2, /3 C = 2, r = —2 and r = 1 as a function of fc in (h 1 .Mpc) 1 
compared to the case with no modification of gravity. 



of gravity is present. When (3 C > 0, Newton's potential shows an anomalous growth with k as soon as small scales 
enter the Compton radius. We also find that the speed of sound is largely affected by /3& implying a shift in the peak 
positions. 

The effect can be seen in the figures 1-5 where three types of couplings have been chosen. When /3 is large and 
/3 C is small, here taken to be /3& = 3 and /3 C = 0, we find that the peaks are shifted due to a change of the sound 
horizon. On the other hand, when (3 C is large and /3& is small, here (3 C = 100 and flt, = 0, the amplitude of the peaks 
is increased. When both are of the same order, the effects are combined together with a new effect coming from 



se = p b p c 



2R $ 
R+ 1 3cf 



(87) 



whose influence can be seen even when /3f, iC are small and no direct deviation on the speed of sound and the amplitude 
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Temperature Fluctuation Temperature Fluctuation 




FIG. 3: The Sachs-Wolfe temperature fluctuation at last scattering for fit, = 3, j3 c = 0, r = —2 and r = 1 as a function of fc in 
(h _1 .Mpc) _1 compared to the case with no modification of gravity. A shift in the position of the peaks due to a modification 
of the sound speed is noticeable. 



Temperature Fluctuation Temperature Fluctuation 




FIG. 4: The Sachs- Wolfe temperature fluctuation at last scattering for /3f, = 0, fi a = 100, r = —2 and r — 1 as a function of 
k in (h _1 .Mpc)~ 1 compared to the case with no modification of gravity. 



of Newton's potential is noticeable. In figure 2 and 5 with j3 c ~p o ~ 2, we show how little effect modified gravity has 
on Newton's potential and the large effect on the CMB peaks from (87). 

In all cases we have chosen m c — mj to enhance the modification of the CMB spectrum. In this case we must 
impose that r < 1.2 and we have chosen two relevant cases with r = 1 and r = —2 where am(a) is either increasing or 
decreasing respectively. When am(a) increases, scales leave the Compton radius eventually. When am(a) decreases, 
scales enter the Compton radius. In the first case, for large enough k, scales are in the Compton radius at horizon 
entry and will eventually leave it. In the latter, scales which are not within the Compton radius at horizon entry will 
eventually enter the Compton radius. This results in different behaviours for Newton's potential as a function of k. 
When /3 C is large enough and am(a) increases, scales with large k tend to increase with a modified rate compared 
to General Relativity for a longer period of time. The charactistic scale k c is k c = 0.006ft, • Mpc -1 when r = 1 and 
k c = O.lh ■ Mpc - when r = —2. In figure 1, one can see the power law growth between k\ and k c when r = —2, 
respectively k c and ki for r = 1, followed by a logarithmic rise for larger values of k. 
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FIG. 5: The Sachs- Wolfe temperature fluctuation Q at last scattering for /3(, = 2, /3 C = 2, r = —2 and r = 1 as a function of fc 
in (h _1 .Mpc) _1 compared to the case with no modification of gravity. 



These examples have been chosen to accentuate the effects. The couplings /3& iC have been taken to be voluntarily too 
large. We expect that the tight bounds from the Planck satellite would lead to stringent constraints on the models. 



C. Engineering modified gravity 

In this section we will show that one can always engineer any modified gravity model characterised by the couplings 
to matter (3b. c and a mass function m(a). This can be achieved using chameleon theories in the cosmological regime 
which arc determined by the couplings ftb,c and the time dependent mass m(a) as a function of the scale factor, as 
long as m > ff implying that the minimum of the effective potential V e f{ is a dynamical attractor. Indeed, the mass 
at the minimum of the effective potential leads to the constraints 



V 



d 2 V 



'PI ~T~ TO P1 

where the couplings /3j can be field dependent. Using the minimum equation 

dV v^„ . Pi 



= -E&^ 



TOpi 



we deduce that the field evolves according to 



U(p _ 3H ^ a A 
~J7 ~ 2 Z,Pi A i ' 



dt m? 

To simplify the analysis, let us assume that Ai = exp(K4/?i0) and /«4^ < 1, we obtain that 

1 



(a) = 



-p(a)da + 4>q 



mpi 7a ini am 2 (a) 

where the sum is over non-relativistic species and we have defined an effective coupling constant /3 as 

fip = E P iPi 

and P = Pi- The minimum equation implies that 



V = V - 3/3 2 



a . . am 2 (a) TOpj 



da. 



(88) 



(89) 



(90) 



(91) 



(92) 



(93) 
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This defines the potential parametrically. Of course, one should check that local tests of gravity are then satisfied with 
this potential. 



Let us give an explicit example: 



D. Power law potentials 



m(a) = m ini (—) r . (94) 

&ini 



In the matter dominated era, we have H(a) = i?mi(^f 1 ) 3 ^ 2 and mini ^> i^ini- We then obtain 



(a) = 4> + 3 f Ani — L (l - (gH»)3+2r) , (95) 
mf-mpi 3 + 2r V v a ' 



As long as 3 + 2r < 0, we can choose <fio such that 



3f3p 

ini J- /-^ini\3+2r 



m? ni m P i 1 3 + 2r| a 
This is a power law behaviour. Similarly we obtain 



(— y +zr (96) 



- a 4 , 3/3 2 p? ni 1 , a mi 
^ ~ A ° + 7< i m| I 6 + 2r ( a } (9?) 



where Ao is the vacuum energy. This leads to a usual chameleon potential where we can identify 

6 + 2r 



3 + 2r 



(98) 



and 



m? ni m P1 6 + 2r TO? ni m P i |3 + 2r| j 



where — 3 < r < —3/2. This defines a chameleonic model coupled to both baryons and CDM particles. When r is 
greater than —3/2, the models are also of the power law type with positive powers of the field. 

Scales such as m ; ~J close to the CMB peak scales and the choice of a power law evolution of the mass function 
(94) between the last scattering surface and now may lead to an incompatibility of the models with both the recent 
growth of large scale structures and/or local gravity tests in the solar system and the laboratory [22]. Models where 
effects of modified gravity could be present both in the recent past of the Universe and at the last scattering surface 
are more likely to require engineering the function m(a) around a ~ 10 and a ~ 1 with different functional forms, 
for instance different power law functions, and interpolating between these two regimes where constraints are scarce. 
This is left for future work 



E. Discussion 



We have shown that for large enough couplings to the baryons and/or CDM, effects on the CMB peaks can be 
induced provided the mass of the scalar field mediating modified gravity is tuned to be close to the scales of the CMB 
peaks. We also know that strong constraints on modified gravity can be obtained at low redshift using large scale 
structures of the Universe [23] . The scales corresponding to the CMB and LSS are different and can be connected by 
an appropriate choice of the mass function m(a). Hence, the CMB and LSS constraints can be seen to be independent 
as they rely on two different ranges of the scale factor for the function m(a). Consequently, unless one can infer a 
function 771(a) from an underlying theory, one must take the constraints given by the position of the CMB peaks and 
the LSS constraints on the growth of structures as independent snapshots of an unknown function 771(a). In fact, at 
the linear level, modified gravity models are characterised by their couplings to matter and the mass function 771(a). 

In this paper we have focused on the CMB physics and its modification due to a scalar field with different couplings 
to baryons and CDM. We have seen that the effect of the baryon coupling is prominent on the speed of sound. In 
particular, this coupling cannot be too large when the mass of the scalar field is of the order of the scale of the CMB 



15 



peaks. For very large couplings, the speed of sound becomes imaginary. For smaller values of the baryon coupling, 
the position of the CMB peaks is shifted due to the alteration of the speed of sound. The CMB coupling influences 
the amplitude of the Newtonian potential. The effect is more relevant when the mass scale am(a) increases as short 
scales enter the Compton radius early and have more time to grow anomalously. 

The subsequent evolution of the m{a) function is not constrained by the CMB physics. It is plausible that 21 cm 
physics may give some indication on the shape of the function m(a) between the CMB and LSS epochs. Constraints 
on structure formation at low redshifts exist already while future surveys will give tighter results. On the whole, the 
landscape of possible modified gravity models, characterised by the couplings to matter and the mass function, is 
almost completely unchartered. In this paper, we have suggested that CMB physics may help unraveling some part 
of its geography. 

IV. CONCLUSION 

Modified gravity models involving a scalar field and generating the acceleration of the Universe seem to require a 
screening mechanism to evade local gravitational tests. The same models must lead to a very slow evolution of the 
scalar field since BBN in order to comply with the stringent bounds on the formation of the elements. In general, 
this results in the fact that these models behave like a A-CDM model since BBN at the background level. At the 
perturbative level and on large scales, the models are characterised by their couplings to matter and the mass function 
m(a), i.e. the mass of the scalar field as a function of time. We have shown that given these inputs one can engineer 
a modified gravity model of the chameleon type (at least when the couplings are constant). We have focussed on 
the CMB peaks and the effect of modified gravity on their structure. We have found that three possible effects can 
occur. For large enough couplings to CDM and a mass scale of order of the sound horizon at last scattering, we have 
a strong modification of Newton's potential and an increase of the CMB peak amplitudes. Far large enough couplings 
to baryons, the speed of sound is reduced implying a shift to large values of the momenta for the CMB peaks. Finally, 
even for moderate couplings, the Sachs- Wolfe temperature could be altered by the combined effect of the couplings 
to baryons and CDM. 

We hope that these effects could be attainable with the forthcoming Planck experiment. Unfortunately, this would 
uncover only one corner of the modified gravity puzzle as large scale structures are sensitive to the mass function at 
much later times. Maybe 21 cm physics could help connecting information from both epochs. This is left for future 
work. 
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